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Abstract — We  have  previously  developed  a  distributed  routing 
algorithm  which  minimizes  the  number  of  overhead  messages  nec¬ 
essary  to  maintain  the  minimum-power  multi-hop  routes  iu  a  mo¬ 
bile  ad-hoc  network,  assuming  a  piece-wise  linear  model  for  the 
motion  of  the  nodes.  In  this  paper  we  extend  the  routing  algo¬ 
rithm  to  Include  clustering  as  well,  to  reduce  further  the  number  of 
overhead  messages  at  the  expense  of  sub-optimal  routes.  A  single 
parameter  controls  the  degree  of  clustering,  and  consequently  the 
degree  of  sub-optimality,  rather  than  arbitrary  parameters  such  as 
maximum  cluster  size  or  maximum  distance  between  nodes.  The 
proposed  algorithm  converges  in  a  finite  number  of  iterations  to 
both  stable  routes  and  stable  clusters,  and  by  setting  the  cluster 
parameter  to  0  collapses  to  the  original  routing  algorithm  with  no 
clusters  and  optimum  routes. 


distributed  algorithm  minimizes  a  global  multi-objective  func¬ 
tion,  which  allows  control  of  the  degree  of  clustering  based  on 
the  minimum-power  routing,  promoting  intelligent  clustering 
by  grouping  into  clusters  nodes  with  similar  trajectories,  rather 
than  arbitrary  parameters  such  as  maximum  cluster  size  or  max¬ 
imum  distance  between  nodes. 

II.  Preliminaries 

This  section  reviews  the  kinetic  routing  algorithm  developed 
in  our  previous  work.  Assuming  a  piece-wise  linear  model  for 
the  motion  of  a  node 


I.  Introduction 

In  a  mobile  ad-hoc  network  (MANET),  many  nodes  are  either 
battery  operated  or  otherwise  power-constrained.  Therefore  it 
is  often  necessary  to  maintain  minimum-power  routes  between 
the  nodes  of  the  network  as  they  move  about,  without  a  signif¬ 
icant  decrease  in  throughput  attributed  to  overhead  messages. 
By  modeling  node  positions  as  a  linear  trajectory,  as  opposed  to 
a  £xed  position  for  a  single  time  instant,  we  have  developed  in 
previous  work  [1]  a  distributed  algorithm  that  renders  a  priori 
(i.e.  at  the  time  of  organization)  the  kinetic  spanning  tree  [2],  or 
the  sequence  of  shortest-path  trees  between  nodes  for  all  time. 
The  algorithm  sustains  virtually  no  increase  in  overhead  mes¬ 
sages  compared  to  organizing  the  routes  for  a  single  time  instant 
through  the  Bellman-Ford  algorithm  [3].  Moreover,  the  kinetic 
tree  requires  minimum  updating  only  when  a  node  changes 
course  since  most  of  the  routes  remain  valid,  rather  than  up¬ 
dating  at  uniform  periods  and  re-propagating  routes  across  the 
whole  network. 

While  the  kinetic  routing  algorithm  minimizes  the  number 
of  overhead  messages,  this  number  may  still  prove  prohibitive 
as  the  size  of  the  network  grows.  A  popular  approach  to  deal 
with  issues  of  scalability  suggests  grouping  nodes  into  clusters 
[4] -[6]:  the  ordinary  nodes  of  a  cluster  elect  a  clusterhead  and 
route  through  it  to  other  nodes  in  the  network,  at  the  expense  of 
sub-optimal  routes.  In  this  paper,  we  propose  the  extension  of 
the  kinetic  routing  algorithm  to  clustering  as  well.  We  study  the 
tradeoff  between  minimum-power  routing  and  clustering,  and 
accordingly  develop  a  kinetic  algorithm  which  simultaneously 
organizes  both  the  clusters  and  the  routes  of  the  network.  The 


Xi{to)  +  Xt{to)-t 
yi{to)  +  yi{to)-t  \  ’ 


(1) 


where  vector  denotes  the  position  of  node  i  at 

the  initial  time  to,  and  vector  2)i(G)]^  its  velocity,  the 

algorithm  computes  the  kinetic  spanning  tree. 

The  squared  distance  between  nodes  i  and  j  is 

Di{t)  = 

—  ttij  t  hij  t  “t“  Cij ,  (2) 


and  the  power  cost  as  a  function  of  time,  required  to  transmit  be¬ 
tween  nodes  i  and  j,  is  de£ned  as  Pij{t)  =  Pji{t)  = 
for  some  constant  n.  The  kinetic  spanning  tree  translates  to  a 
forwarding  table  at  each  node  whose  entries  are  the  minimum- 
power  paths  in  time  to  the  sink  node.  The  table  of  node  i  ap¬ 
pears  in  Fig.  1.  The  shaded  area  indicates  its  composite  power 
cost  P,{t)  =  mm{Py^]s(t),  required  to 

route  from  i  to  S'  through  the  next-hop  nodes  ji,  j2.  and  j^',  [] 
denotes  the  ordered  nodes  comprising  the  rest  of  the  path.  Node 
i  forwards  to  node  ji  for  0  <  f  <  G ,  to  node  j2  for  G  <  t  <  t2, 
to  node  jo  for  t2  <  t  <  to,  and  to  node  j2  for  to  <  t. 

As  a  simple  example,  Table  1(a)  lists  the  trajectories  of  four 
nodes  routing  directly  to  the  sink  node  S  at  time  to  =  0,  except 
D  A.  The  changes  in  the  kinetic  tree,  computed  a  priori  at 
time  to  =  0,  appear  in  Table  1(b). 


III.  Interaction  of  Routing  and  Clustering 
A.  Power  Function 


OThis  work  was  supported  by  the  Advanced  Research  and  Development  Ac-  The  kinetic  routing  algorithm  described  in  the  previous  sec- 

tivity  (ARDA)  under  contract  number  706400.  tion  effectively  minimizes  at  node  i  the  power  function 
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Fig.  1.  New  and  current  power  costs  in  the  pruning  step  at  node  j. 

TABLE  I 

Changes  in  the  kinetic  spanning  tree. 


The  values  of  (m)  and  Rm  {m)  depend  on  the  chosen  clus¬ 
tering  protocol.  Here  we  assume  that  ordinary  nodes  route  di¬ 
rectly  (i.e.  single-hop)  to  their  clusterhead,  while  the  cluster- 
heads  route  through  each  other  in  a  multi-hop  fashion  to  the 
sink  node.  Also  the  sink  node  will  always  choose  to  stay  a  clus¬ 
terhead  to  minimize  the  power  function. 

When  an  ordinary  node  changes  course,  the  features  used  to 
group  it  with  its  cluster  change,  and  so  it  may  no  longer  belong 
to  that  cluster;  accordingly,  we  reinitialize  this  node  as  a  single 
clusterhead  and  allow  it  to  reorganize.  If  a  clusterhead  changes 
course,  we  reinitialize  it  as  a  single  clusterhead;  in  addition, 
the  ordinary  nodes  associated  with  this  cluster  no  longer  have  a 
route  to  the  node,  and  so  must  reinitialize  as  single  clusterheads 
and  be  allowed  to  reorganize. 

The  message  complexity  incurred  to  reorganize  the  multi¬ 
hop  routes  of  n  clusterheads  to  the  sink  varies  as  0{n?)  [7], 
hence  for  this  particular  protocol  Rn{m)  =  (n  -f  m  —  1)^  and 
Rm{m)  =  (n-fl)^.  Substituting  n  =  N/m  into  (4)  and  solving 
yields  the  solution  for  m*  which  minimizes  R{m) 


t  Change 
T252  A-+B 
2.193  D^C 
3.754  A  ->  S 

4.085  B  ^  A 

5.692  C  ^  A 

8.378  D  ^  A 

8.452  D  ^  S 

(b) 


OO 

Fr  =  J  m)  dt  (3) 

^0 

in  a  distributed  manner,  where  to  denotes  the  execution  time. 

If  clusterhead  i  joins  another  cluster  j  as  an  ordinary  node, 
node  i  then  routes  to  S  through  node  j.  Accordingly  its  com¬ 
posite  cost  Pi{t)  increases  to  P/(f)  =  Pij{t)  +  Pj{t)  forcing  an 
increase  in  Ff  as  well,  hence  sacrifcing  the  minimum-power 
routes  to  reduce  the  overhead  messages  necessary  to  reorganize 
the  network  in  the  event  of  a  trajectory  change.  Note  that  if  node 
j  was  routing  through  i  before  the  clustering,  that  path  no  longer 
exists  after  the  clustering,  hence  Pj^  (as  well  as  the  power  func¬ 
tion  of  other  clusterheads  routing  through  i)  increases  as  well. 

B.  Reorganization  Function 

Given  N  nodes  in  a  network  and  assuming  on  average  m 
nodes  per  cluster,  we  determine  the  optimal  value  m*  which 
minimizes  the  cost  of  reorganizing  the  network  in  the  event 
of  a  course  change.  Let  n  indicate  the  number  of  clusters, 
then  pn  =  n/N  represents  the  probability  that  a  clusterhead 
changed  trajectory  conditioned  on  some  node  changing  trajec¬ 
tory,  with  Rn{rn)  the  associated  cost  of  reorganization.  Like¬ 
wise  Pm  =  {N  —  n)  /N  represents  the  conditional  probability 
that  an  ordinary  node  changed  trajectory,  with  associated  cost 
Rm{m).  We  formulate  the  conditonal  cost  of  reorganization  as 

R{m)  =  Pn-Rnim) +Pm-R7n{rn).  (4) 


.  N(N  -  4) 

fYl  ^  _ 

X  =  Aft  [27 N'^  -  2A8N  -I-  432  -I-  3(Ar  -  4)  yJ‘i{27N'^  -  2mN  +  432)] 

The  reorganization  function  below  quanti£es  the  degree  to 
which  node  i  departs  from  the  optimal  con£guration:  R{m)  is 
much  steeper  farther  from  m*  and  so  a  smaller  cluster  will  ben- 
e£t  more  from  accepting  an  additional  node  than  a  larger  cluster. 

In  minimizing  the  function  to  0,  the  network  converges  to  m* 
nodes  per  cluster. 

max{mj  ,m*  } 

Ff^  =  J  R{m)  dm  (6) 

C.  Multi-Objective  Function 

Before  organization,  each  node  in  the  network  begins  as  a 
single  clusterhead.  If  the  network  organizes  through  minimiz¬ 
ing  the  power  function  alone,  no  clustering  takes  place  and  so 
constructs  the  kinetic  minimum-power  routes  between  the  clus¬ 
terheads  to  the  sink.  At  the  other  extreme,  minimizing  only  the 
reorganziation  function  results  in  a  network  with  n*  =  N/m* 
clusters  to  reduce  the  number  of  overhead  messages  to  maintain 
the  routes.  In  general,  we  seek  a  compromise  between  the  two 
opposing  functions  by  minimizing  the  multi-objective  function 

F,  =  F[  +  5-F^,  (7) 

where  the  only  parameter  5  controls  the  degree  of  clustering 
suited  to  a  particular  application,  rather  than  arbitrary  param¬ 
eters  such  as  maximum  cluster  size  or  maximum  distance  be¬ 
tween  nodes  [4]-[6].  F/^  forces  the  merging  of  nodes,  and  sub¬ 
ject  to  this  constraint  F//^  gives  preference  to  clustering  those 
nodes  moving  with  similar  trajectories,  since  the  power  cost  in 
time  between  them  will  be  smallest.  So  the  power  function  will 
cluster  in  an  intelligent  manner  based  on  the  feature  of  similar 
trajectories.  We  can  imagine  separate  troops  of  soldiers  moving 
in  a  military  operation,  all  soldiers  in  a  troop  following  similar 


Node 

Trajectory 

S 

(0,0) 

A 

(-1  -O.lf,  1) 

B 

(-0.4t,  1.5-0.2t) 

C 

(2  -  0.5L  1 -F  0.2i) 

D 

(0,2-F0.3t) 

(a) 
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trajectories.  Our  proposed  clustering  based  on  minimum-power 
routing  will  favor  making  each  troop  a  separate  cluster,  espe¬ 
cially  if  there  is  a  great  distance  between  troops.  Therefore  the 
clustering  can  reveal  the  structure  behind  the  network,  if  one 
exists. 

As  pointed  out  in  Section  III- A,  when  a  node  i  joins  another 
cluster,  it  alters  not  only  the  local  objective  function  Fi,  but 
also  the  function  of  other  nodes  in  the  network.  Hence  we  min¬ 
imize  the  global  multi-objective  function  F  =  Fi  through 
a  distributed  technique  known  as  “region  competition”  [8].  The 
following  section  provides  the  details  of  the  algorithm. 

IV.  Kinetic  Routing  and  Clustering  Algorithm 

Initialization 

Each  node  maintains  a  forwarding  table  for  itself  and  its  neigh¬ 
boring  clusterheads  (not  ordinary  nodes).  It  also  maintains  the 
number  of  nodes  in  each  of  their  clusters.  Each  node  j  in 
the  network  begins  as  a  clusterhead.  It  adds  to  its  table  j  the 
power  cost  to  the  sink,  Pjs{t).  If  S  is  not  a  neighbor  of  j,  then 
Pjs{t)  =  oo.  Node  j  distributes  Pjs{t)  to  its  neighbors. 
Iteration  Step 

1.  UPDATE:  At  the  current  time  step,  node  i  receives  Pj^]g{t) 
from  neighboring  clusterheads  j  €  ji,  •  •  • ,  and  adds  Pj\js  to 
its  table  j. 

2.  PRUNE  (clusterheads  only):  Node  i  computes  the  current 
cost  Pij[]g{t)  =  Pij{t)  +  Pj[jg{t)  for  each  node  j.  The  cur¬ 
rent  costs  and  the  previous  costs  (i.e.  those  in  table  i  from  the 
previous  step)  are  compared  amongst  each  other.  Only  those 
(minimum)  costs  which  contribute  to  the  composite  cost  Pi{t) 
are  retained  in  table  i.  Eig.  1  displays  both  the  current  and  pre¬ 
vious  costs  in  table  i  at  the  current  step.  They  appear  as  six 
parabolic  functions  in  time,  however  only  three  of  them,  namely 
Pih[]SiPih[]S,  and  Pij^[]g,  contribute  to  P,(t).  Only  the  cur¬ 
rent  costs  added  to  the  table  in  the  current  step  are  then  dis¬ 
tributed  to  neighboring  nodes,  since  the  previous  costs  added  to 
the  table  were  already  transmitted  in  those  steps,  and  need  not 
be  retransmitted. 

3.  CLUSTER  (single  clusterheads  and  ordinary  nodes  only): 
At  this  step  node  i  may  join  a  neighboring  cluster  j,  if  such 
an  action  lowers  the  global  multi-objective  function  of  the  net¬ 
work,  which  translates  to  a  negative  change  over  the  local  multi¬ 
objective  functions  of  those  nodes  affected  by  the  clustering.  In 
region  competition,  node  i  computes  for  each  of  its  n  neigh¬ 
boring  clusterheads  j  G  ji,  j2,  ■  ■  ■ ,  jn  the  change  in  the  global 
function  AF[j]  which  results  from  joining  cluster  j.  i  then  joins 
the  cluster  j*  which  offers  the  most  negative  change.  If  there  is 
no  negative  change  amongst  the  candidate  neighboring  cluster- 
heads,  node  i  remains  in  its  current  state. 

If  the  clustering  takes  place,  i  must  be  removed  from  the  ta¬ 
bles  of  all  neighboring  clusterheads  which  route  through  it  to  the 
sink,  resulting  in  the  modi£ed  Pj_^  (t) ,  (^)  >  •  •  • )  Pj^  {t)  which 

can  be  used  to  compute  the  changes  in  the  local  functions  of 
these  affected  nodes.  Since  i  contains  the  routing  table  of  all  its 
neighboring  clusterheads,  i  can  compute  these  changes  locally 
without  any  additional  exchange  of  messages^.  The  change  in 

^  An  exact  calculation  would  involve  removing  i  from  the  tables  of  all  clus¬ 
terheads  in  the  network,  not  just  the  neighboring,  but  this  would  not  allow  local 
computation  of  the  change;  moreover,  the  neighboring  nodes  are  those  most 
affected  by  the  clustering. 


the  local  function  of  tcalculation  affected  i  is  computed  using 
the  composite  cost  P/(f)  =  Pij{t)  +  Pj{t)  from  routing  di¬ 
rectly  through  clusterhead  j.  In  addition,  one  needs  to  consider 
the  new  m'  =  m'  =  mj  +  1  and  =  to'  of  the  ordinary 
nodes  of  cluster  j,  as  well  as  the  corresponding  changes  in  the 
number  of  nodes  in  the  cluster  that  i  leaves  if  it  joins  If  i 
joins  j*,  it  does  so  by  sending  one  join  message  to  its  neighbors 
which  achieves  three  objectives:  I.  it  informs  j*  that  it  will  be 
joining  it;  2.  it  updates  the  value  of  to  of  the  ordinary  nodes  in 
cluster  j*  and  the  previous  cluster  of  i;  3.  it  removes  any  entries 
on  the  tables  its  neighbors  maintain  that  contain  i  on  its  path. 

The  reader  should  be  aware  of  the  special  cases  in  which  a 
join  message  is  denied: 

1.  Node  i  attempts  to  join  clusterhead  j*,  which  itself  has  joined 
another  clusterhead  as  an  ordinary  node  at  the  same  iteration. 
In  this  case  j*  sends  a  deny  message  and  node  i  remains  in  its 
current  state. 

2.  Two  single  clusterheads  wish  to  join  each  other’s  cluster.  In 
this  case  the  node  with  the  smaller  ID  (or  some  other  convention 
may  be  used)  joins  the  cluster  of  the  node  with  the  larger  ID.  No 
deny  message  is  required. 

Given  the  non-convexity  of  the  global  multi-objective  func¬ 
tion  for  (5  >  0,  simulated  annealing  is  introduced  to  avoid  con¬ 
vergence  to  a  local  minimum.  Note  that  with  (5  =  0,  the  kinetic 
algorithm  simply  collapses  to  the  kinetic  routing  algorithm. 

V.  Modified  Power  Cost 

Employing  the  power  cost  described  thus  far  clusters  based 
on  the  state  of  the  network  as  f  ^  oo,  since  the  last  entry  in 
the  forwarding  table  of  a  node  approaches  oo  as  f  ^  oo  (see 
Pjj2[]s(f)  of  Eig.  1).  This  attenuates  the  contribution  of  the 
other  entries  to  Pi{t)  and  to  Ff .  To  avoid  this  problem,  we 
bias  the  clustering  towards  to  since  as  f  ^  oo  all  the  nodes  will 
have  changed  trajectory  anyway.  We  de£ne 

K(f)  =  (8) 

as  the  probability  that  a  node  i  is  continuing  on  its  present  tra¬ 
jectory,  where  the  Poisson  parameter  indicates  the  average 
time  the  node  follows  a  course,  and  U  the  time  its  current  tra¬ 
jectory  began. 

Assuming  independent  node  trajectories,  Pij{t)  =  Pi{t)  ■ 
Pj  (f)  describes  the  probability  that  nodes  i  and  j  are  continuing 
on  their  respective  courses  at  time  t.  The  modi£ed  power  cost 
below  probabilistically  weights  Pij  (f)  to  reCtect  the  uncertainty 
of  the  link’s  existence. 


Pij  (t) 


Pijjt) 
Pij  it) 


(9) 


A  low  modi£ed  power  cost  favors  a  low  power  cost  with  high 
probability.  It  readily  substitutes  in  the  algorithm  described  in 


^For  a  correct  implementation  of  region  competition,  the  two  components  of 
the  change  in  the  global  function  AF^]  must  be  normalized  in  the  sum  by  the 

largest  respective  components  of  the  candidates:  — ^  [j]  | } 

max  {|AF-^^[j]|}  ’ 
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the  previous  section  by  simply  changing  the  alegbraic  sum  to 
the  sequel: 


PijQsit) 


Pij{t)  0  Pj[]s{t)  — 


p. 

-*7,1 


tj  (t) 


Pijt)  ■Pjl]s{t) 
■PjQsit) 


(10) 


[1S0[]S 


Pi+P 


J[1S 


{aij  +  aj[]s)P  +  {bij  +  bj[]s)t  +  {cij  +  Cj[]5) 
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Now  the  £ve  parameters  Cij[]s,  Pij[]S,  and 

describe  Pij^s{t)  when  sending  a  routing  message,  rather  than 
three  for  P^j[]s{t). 

It  can  be  shown  that  that  and  in  Section  III-B  are 
still  valid  under  the  assumption  of  Poisson-distributed  trajectory 
changes,  and  so  the  associated  equations  presented  to  generate 
the  reorganization  function  still  prove  valid. 

VI.  Simulation  Results 

The  following  section  contains  simulation  results  for  the  rout¬ 
ing  and  clustering  algorithm  for  a  network  of  N  nodes  uni¬ 
formly  spread  over  a  10  x  10  mile  area  at  initialization  time 
t  =  0.  N  varies  as  {20,40,80}  by  simulation,  and  all  results 
were  averaged  over  10  trials.  We  assign  each  node  a  random  ve¬ 
locity,  whose  speed  does  not  exceed  60  miles/hour.  The  trans¬ 
mission  radius  of  the  each  node  is  limited  to  4  miles.  We  as¬ 
sume  that  the  underlying  communication  links  requires  100  ms 
to  transmit  a  packet  across  single  link,  and  set  the  unit  time  for 
one  iteration  accordingly.  The  course  changes  of  the  nodes  are 
Poisson  distributed  with  value  (3  =  0.03,  hence  a  node  changes 
its  trajectory  on  average  every  33s. 


Fig.  2.  Number  of  clusters  formed  vs.  the  clustering  parameter  with  £xed 
trajectories. 


Fig.  3.  Number  of  routing  messages  vs.  the  clustering  parameter  with  £xed 
trajectories. 


A.  Fixed  Trajectories 

The  experiments  with  £xed  trajectories  allow  the  reader  to 
extract  the  intrinsic  properties  of  the  network  characterized  by 
the  clustering  parameter,  namely  the  number  of  clusters  formed, 
the  number  of  routing  messages  required  for  self-organization, 
and  the  transmission  power  of  the  associated  routes.  Each  sim¬ 
ulation  was  run  for  30  seconds,  terminating  before  any  expected 
trajectory  changes. 

An  increase  of  the  clustering  parameter  from  the  value  0 
forces  the  reorganization  function  to  play  a  more  signi£cant  role 
in  the  multi-objective  function,  emphasizing  the  cost  of  routing 
messages  to  organize  the  network.  This  attracts  clusterheads  to 
join  other  clusters  as  ordinary  nodes,  reducing  the  total  number 
of  clusterheads  and  so  the  0{rP)  complexity  of  the  algorithm 
accordingly.  Fig.  2  illustrates  the  structure  of  the  organized  net¬ 
work  for  the  values  5  G  (0.0, 0.2,  0.5, 0.7, 1.0}.  The  optimal 
number  of  clusterheads  for  N  G  {20, 40, 80}  given  through  (5) 
is  n*  G  {1.72,  2.23, 2.90}  respectively,  and  we  see  in  each  case 
that  the  number  of  clusters  approaches  n*  as  S  increases. 

For  the  same  operation  points  of  the  clustering  parameter. 
Fig.  3  displays  the  corresponding  number  of  routing  messages 
required  for  self-organization  of  the  network  from  initialization 
to  convergence  to  stable  routes.  The  savings  in  the  number  of 
messages  with  increasing  6  is  more  pronounced  as  the  size  of 


the  network  grows,  from  a  factor  of  23%  with  iV  =  20  to  a 
factor  of  448%  with  N  =  80.  It  is  worth  mentioning  that  the 
iteration  complexity  of  the  algorithm,  hence  the  number  of  iter¬ 
ations  required  for  convergence,  varies  in  proportion  to  the  mes¬ 
sage  complexity.  For  <5  =  0  and  N  G  {20, 40, 80}  the  algorithm 
requires  on  average  {5, 10,  20}  iterations  for  convergence,  each 
decreasing  with  larger  values  of  the  clustering  parameter. 

Fig.  4  displays  the  sum  of  the  power  function  over  the  N 
nodes  in  the  network  versus  <5,  normalized  by  this  quantity  for 
(5  =  0;  hence  the  transmission  power  associated  with  the  min¬ 
imum  routes  assumes  unity.  The  plot  illustrates  the  departure 
from  the  minimum-power  routing  as  5  increases  from  0.  As  ex¬ 
pected,  the  power  increase  for  routing  as  <5  increases  becomes 
more  pronounced  for  larger  values  of  JV.  For  <5=1,  the  increase 
is  only  55%  for  N  =  20,  but  grows  to  545%  for  N  =  80. 

These  results  show  the  tradeoff  between  minimum-power  and 
the  number  of  routing  messages,  suggesting  the  operation  of 
the  network  at  an  application-suited  point  characterized  by  the 
value  of  5. 

B.  Randomly  Changing  Trajectories 

The  simulations  presented  in  this  subsection  are  identical  to 
those  in  the  previous  subsection,  except  that  here  we  run  them 
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Fig.  4.  Normalized  power  function  vs.  the  clustering  parameter  with  £xed 
trajectories. 


Fig.  5.  Number  of  routing  messages  vs.  the  clustering  parameter  with  changing 
trajectories. 

for  an  extended  dnration  of  £ve  minutes,  and  moreover  the 
nodes  are  allowed  to  change  trajectory  according  to  the  Poisson 
distribntion  with  parameter  f3  =  0.03,  and  reorganize  accord¬ 
ingly.  We  present  the  same  format  of  results  in  the  nnmber  of 
messages  required  to  organize  and  maintain  the  routes  through¬ 
out  the  simulation  period,  and  the  associated  normalized  power 
function  over  all  the  nodes  in  the  network.  The  normalized 
power  function  is  summed  from  the  intialization  time  t  =  0 
to  t  =  300s.  We  cannot  report  the  number  of  clusters  formed 
throughout  the  run  since  this  value  ttuctuates  with  every  course 
change. 

Fig.  5  displays  an  even  more  dramatic  reduction  in  the  num¬ 
ber  of  routing  messages  required  in  the  presence  of  trajectory 
changes  as  6  increases.  Unfortunately  such  a  reduction  results 
in  a  corresponding  drastic  increase  in  the  power  requirements 
for  transmitting  from  all  nodes  in  the  network  to  the  sink  with 
sub-optimal  routes,  as  reClected  in  Fig.  6. 

VII.  Conclusions  and  Current  Work 

We  have  described  a  distributed  algorithm  which  simultane¬ 
ously  organizes  the  routing  and  clustering  of  a  network,  exploit¬ 
ing  the  theory  of  kinetic  spanning  trees.  As  in  kinetic  routing 
algorithm  introduced  previously,  the  proposed  algorithm  con¬ 
verges  in  a  £nite  number  of  iterations  to  both  stable  clusters  and 


Fig.  6.  Normalized  power  function  vs.  the  clustering  parameter  with  changing 
trajectories. 

stable  multi-hop  routes  between  the  clusterheads.  Moreover,  the 
organization  updates  are  event-driven  by  a  change  in  the  trajec¬ 
tory  of  a  node,  as  opposed  to  updating  at  uniform  periods  as  in 
conventional  algorithms. 

The  simulation  results  show  that  by  increasing  the  clustering 
factor  6,  we  can  signi£cantly  decrease  the  number  of  messages 
necessary  to  organize  and  maintain  the  routes  from  all  the  nodes 
in  the  network  to  the  sink  node.  This  reduction,  however,  comes 
at  the  expense  of  sub-optimal  routes  in  terms  of  power  tranmis- 
sion. 

Our  current  work  involves  simulating  networks  where  the 
nodes  move  about  as  groups  according  to  similar  trajectories, 
rather  than  independently  according  to  some  random  distribu¬ 
tion.  The  simulations  should  illustrate  the  true  utility  of  orga¬ 
nizing  the  network  by  minimizing  the  multi-objective  function, 
where  the  groups  are  formed  intelligently  into  functional  clus¬ 
ters. 
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